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Abstract. Proving the existence of a solution to a system of real equations is 
a central issue in numerical analysis. In many situations, the system of equa- 
tions depend on parameters which are not exactly known. It is then natural 
to aim proving the existence of a solution for all values of these parameters in 
some given domains. This is the aim of the parametrization of existence tests. 
A new parametric existence test based on the Hansen-Scngupta operator is 
presented and compared to a similar one based on the Krawczyk operator. It 
is used as a basis of a fixed point iteration dedicated to rigorous sensibility 
analysis of parametric systems of equations. 



Notations 

Vectors are denoted by boldface symbols, and interval, interval vectors and ma- 
trices by bracketed symbols. Let f : E — > F be a function, and X C dom(f). Then 
f (X) := {f (x) e F|x e X} is the range of f on X. 

1. The Hansen-Sengupta Existence Test 

The presentation given here follows the one given by Neumaier in [4]. The 
interval Gauss-Seidel is defined as follows: First in dimension one, 

(1) [7] ([a], [6], [x]) := a{x £ [x] \ 3a £ [a],3b G [b],ax = b}. 

In the case where ^ [a], one obtains the expression [7] ([a], [b], [x]) = ([&]/[«]) H [x] 
(cf. [4] for the expression in the case G [a]). Then, the multidimensional Gauss- 
Seidel is then defined as follows: [r]([j4], [b], [x], [z]) [x'] where 

(2) K]:=[7](K], , n). 

j<i j>i 

Remark 1. In the traditional definition of the interval Gauss-Seidel operator, the 
interval vector [z] is equal to [x] (and hence does not appear explicitly in its def- 
inition). Using instead [z] = R" disactivates the intersection with the previous 
domain and can be useful for some applications (cf. Section [3]). 

Then, the Hansen-Sengupta operator [5] can be defined as follows 

(3) x-f [r]f[X],-[y],[x]-i,[z]-i). 



where [X] e IM"''" and [x], [y], [z] e IE". The following theorem shows how the 
Hansen-Sengupta operator can be used to improve the enclosure and prove the 
existence of solutions (cf. [4]). 
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Theorem 1. Let [x], [y], [z] e IM", x e [x] and [X] e IM"''" such that: [x] C [z], 
f(i) e [y] and [X] D {^(x) £ R"><"|x € [x]}. // [x'] denotes @ i/ien; 

(1) X S [x] anrf f(x) = implies x G [x']. 

(2) If ^ ^ [x'] C int[x] f/ien f /las an unique zero in [x']. 

Remark 2. The interval evaluation of the derivatives can be replaced by Lipschitz 
interval matrices to release the differentiability hypothesis, and to slope matrices 
to improve the enclosure (though uniqueness of the solution is lost when slopes are 
used, cf. [5 for details). 

A preconditioning is usually coupled to this kind of operator: The preconditioned 
system C ■ f(x) = 0, where C S K"^" is nonsingular, is equivalently solved instead 
of f(x) = 0. The preconditioning matrix C is chosen so that C • f is close to the 
identity in the domain considered, hence improving the efficiency of the operator. 

2. Parametrization of the Hansen-Sengupta Existence Test 

Functions with parameters arc considered in this section. Let f : Rp x M" — > R" 
be a function of n variables and p parameters. Parameters will be denoted by the 
vector a and variables by the vector x. The parametric Hanscn-Scngupta operator 
is expressed applying its non-parametric version to different inputs. A more general 
parametric Hansen-Sengupta (which was dedicated to quantified parameters thanks 
to the usage of the Kauchcr arithmetic) was proposed and used in [1] . 

Theorem 2. Let [x], [y], [z] G IR", [a] G IR^, ii: G [x] and [X] G IR"''" such that: 
[x] C [z], f([a],i) C [y] and [X] D {^(a,x) G R"X"ia G [a],x G [x]}. // [x'] 
denotes ([3]) then: 

(1) a G [a] and x G [x] and f(a, x) = implies x G [x']. 

(2) If ^ ^ [x'] C int[x] then for every a G [a], f(a, •) has an unique zero in [x']. 

Proof. Fix an arbitrary a G [a] and define g : R" — * M" by g(x) = f(a, x). 
We are going to apply Theorem [1] to g. First, g(x) = f(a, x) G [y]. Second, as 
g(x) = f(a,x), 

(4) |^(x) G M"X" I X G [x]| = |^(a,x) G R"><" | x G [x]| C [X]. 

Therefore, Theorem [1] can be applied to g and the domain [x], and shows that if 
[x'] denotes ^ then 

(1) g(x) = (that is f(a,x) = 0) implies x G [x']. 

(2) ^ [x'] C int[x] implies the existence of an unique zero of g (that is of 
f(a,.)) in[x']. 

This holds for every a G [a] and hence concludes the proof. □ 

An interval extension of f can be used to compute [y] satisfying f([a],x) C [y]. 
Using the mean- value extension to compute [y] and the usual inverse midpoint pre- 
conditioning gives rise the following parametric Hansen-Sengupta operator, denoted 
by [i7]f_[a]([x], [z]) in the sequel: 

(5) i+[r](c-[X] , -C-[f](&,i)-(C-[A])-([a]-a) , [x]-i), 

with [X] = [f ]([a],[x]), [A] = [£]([a],x), C = {m\A[X])-\ i = mid[x] and 
a = mid [a]. The operator ^ will be denoted by [i?]fja]([x], [z]) in the sequel. 
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Experiments presented in the sequel will be carried out using the natural interval 
extensions of f , ^ and 

Remark 3. In ([5]), the expression C ■ [f](a, x) + (C • [A]) ■ ([a] — a) is used instead 
of C ■ ([f](a, x) + [A] ■ ([a] — a)) because it is a more efficient interval evaluation. 

The preconditioned parametric Hansen-Sengupta operator is compared to the 
preconditioned parametric Krawczyk operator 

(6) [Al^[a]([x]) := 5c+[Al(c.[X],-[y],[x]-ic) 

(7) [Al([A],[b],[x]) [b] + (/-[A]).[x]. 

proposed in [S], where the same interval enclosure [y] := C • [f](a, x) + (C • [Afj ■ 
([a] ~a) of {C-f(a, x)|a £ [a]} is used in ([6]) and in ^ (this point is not detailed in 
i). As in dsn, [X] = [^]([a], [x]), [A] = [£]([a],i), C = (mid[X])-i, x = mid[x] 
and a = mid [a]. Theorem [21 the inclusion [ATjf [a]([x], [z]) implies the existence of 
an unique solution to each system f (a, •) = for a G [a]. 



Example 1. Let 




[a] := ([i-e,^ + e],[-e,e],[l-e,l + e]) and [x] ([-0.2, 0.2], [-0.7, l-l])- The 
uncertainties width is set to e = 0.025. The set { x G [x] | 3a G [a],f(a, x) = } 
is approximately represented on the left hand side graphic of Figure [T] solving the 
2x2 system of equations for a finite set of parameters values inside [a] . 

Both operators ([5]) and ([6| are used to improve the initial enclosure [x] by com- 
puting the sequences 

(9) [x,+i] = [F]f,[,]([x,],[xfe]) 

(10) [y^+i] = [A1f,[a]([yfc]) 

with [xq] = [yo] = [x]. The following table summarizes the results obtained with 
both operators: 





Final enclosure 


Existence step 


Hansen-Sengupta 


([-0.074,0.075], [0.831,0.901]) 


3 


Krawczyk 


([-0.074,0.075], [0.831,0.901]) 


2 



The final enclosure is also shown on the left hand side graphic of Figure [T] These 
results seem to show that the parametric Krawczyk operator is sharper than the 
parametric Hansen-Sengupta operator: they both compute the same final enclosure 
while the former proves the existence one step before. This is surprising since in 
their non parametric form the Hansen-Sengupta operator is proved to be sharper in 
general than the Krawczyk operator (cf. [4]). However, a closer study shows that 
the Hansen-Sengupta operator is actually sharper: The right hand side graphic of 
Figure [1] shows the ratio 

||wid([yfe])||-wid([x,])|| ^ 
wid([x/c])|| 

As this graphic shows, the enclosure computed by the Hansen-Sengupta operator 
is alway sharper. The difference is sensible at the first iterations (reaching approx- 
imately 20% at step 5), and converges to as goes to infinity (the dashed line 
corresponds to 12 exp(— 0.46fc) for information about the convergence rate to 0). 
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Figure 1. Left: Approximate solution set and its bounding box 
computed in Example [T] Right: Ratio of the enclosures widths for 
the parametric Hansen-Sengupta and Krawczyk operators plotted 
for each step. 



This also explains why the existence proof is posponed of one step for the Hansen- 
Sengupta operator: the enclosure computed at step = 1 by this latter operator is 
too sharp to obtain [X2] C int [xi]. 



3. Rigorous Sensitivity Analysis 

A direct application of Theorem [2] requires an initial domain. However, in prac- 
tice this initial domain is often not available. Instead, an approximate solution x* 
for a nominal parameter value a* G [a] is available. In the sequel, x* is supposed 
to satisfy exactly f(a*,x*) = 0, but the usage of an approximate solution has no 
incidence in practice. From the sensitivity analysis point of view, we need to prove 
that each parameter a G [a] is mapped to an unique solution x and to enclose the 
set of these solutions. 

Provided that the parameters domains are small enough, the iteration 

(12) [xfc+i] = [i?]([xfe],R") ; [xo] =x* 

will converge to [xqo] which satisfies [xqo] ~ [if]([xoo], K"). Although this limit can 
be proved to contain the solution set, the inclusion [x^+i] C int[xfc] will never be 
satisfies because this iteration somehow translates and infiates the initial approxi- 
mation X* . It is more practical to additionally inflate each iterate of a fixed ration 
6 so as to obtain the inclusion [x^^+i] C int[xfc] after a finite number of steps, and 
hence prove the inclusion of the solution set. This leads to Algorithm [TJ 

Remark 4. Fixing a maximum number of steps /cmax forces the termination of the 
algorithm. Smarter stopping criteria can easily be implemented. Also, it can be 
noted that once the existence proof has succeeded, the iteration becomes contracting 
and encloses the solution set. Therefore, Line [1] can be replaced by [x'] ^ [x] once 
success is true. 

As shown by the next examples, Algorithm[T]can be used as a rigorous sensibility 
analysis of the solution to a parametric system of equation: Being simply given an 
approximate solution for a nominal parameter value. Algorithm [T] allows rigorously 
bounding the variations of the solutions w.r.t. the variations of parameters. 

Example 2. Let f and [a] be defined as in Example [l] and consider the approximate 
solution X* :— (0.01,0.85), which is represented by a cross in Figure [2l The first 



SENSITIVITY ANALYSIS USING A FIXED POINT INTERVAL ITERATION 



5 



Algorithm 1: 


Input: f : RP X R" ^ M", [a] G IW, x* e K" 




Output: [x] e M" 




1 k 


max ^ 10; /* Maximum number of iterations 


*/ 


2 5 


^ 1.01; /* Inter-step inflation ratio 


*/ 


3 success ^ false; 




4 [x] ^x*; 




5 repeat 




6 


[x'] ^ mid[x] + (5 ([x] - mid[x]); 




7 


[x]^[i/]f.[,]([x'],K"); 




8 


if ( [x] C int[x'] ) then success ^ true; 




9 until ( k > fcniax ) ; 




10 if ( -1 success ) then [x] ^ M"; 




11 return ( [x] ); 
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Figure 2. Approximate solution set and its bounding box com- 
puted in Example [21 The initial approximate solution is repre- 
sented by a cross and the boxes computed by the first steps in 
dashed lines. 

steps of Algorithm [1] are represented by dashed boxes on Figure [2] The existence 
is proved at after four iterations, and hence the solution set is enclosed. 

Related Work. In addition to [5| , we have found two references [3l [6] which use 
interval analysis for rigorous sensitivity analysis. The advantage of the method 
presented here is that it does not require an initial enclosure of the variations. 
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